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Abstract. We present here a study of chemical abundances in the giant elliptical galaxy M 87 using high-resolution spectra 
obtained with the Reflection Grating Spectrometers during two deep XMM-Newton observations. While we confirm the two- 
temperature structure of the inter-stellar medium (ISM) in M 87, we also show that a continuous temperature distribution 
describes the data equally well. The high statistics allows us, for the first time, to determine relatively accurate abundance 
values also for carbon and nitrogen. The comparison of the abundance ratios of C, N, O and Fe in the ISM of M 87 with those 
in the stellar population of our Galaxy shows that the relative contribution of core-collapse supernovae to the enrichment of the 
ISM in M 87 is significantly less than in the Milky Way and indicates that the enrichment of the ISM by iron through Type la 
supernovae and by carbon and nitrogen is occurring in parallel. This suggests that the main source of carbon and nitrogen in 
M 87 are the low- and intermediate-mass asymptotic giant branch stars. From the oxygen to iron abundance ratio in the hot gas 
we estimate that the relative number of core collapse and type la supernovae contributing to the enrichment of the ISM in the 
core of M 87 is ~ 60% and ~ 40% respectively. The spatial distributions of iron and oxygen are different. While the oxygen 
abundance distribution is flat the iron abundance peaks in the core and has a gradient throughout the 4' wide field of view of the 
instrument, suggesting an early enrichment by core-collapse supernovae and a continuous contribution of type la supernovae. 
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1. Introduction 

The giant elliptical galaxy M 87 resides at the center of the 
nearby Virgo Cluster. Its relatively small distance of ~ 16 Mpc 
makes it a perfect target to study processes that are important 
in all clusters of galaxies, but to a much greater detail than pos- 
sible elsewhere. Among such studies, the investigation of the 
structure and physics of cooling cores (previously thought to 
harbour cooling flows) and of the heavy element enrichment of 
the intra-cluster medium (ICM) gain a central importance. This 
is due to the fact that M 87 is on one hand a famous cooling 
core cluster and on the other hand provides with an ICM tem- 
perature range of about 1-3 keV an X-ray spectrum very rich 
in metal emission lines that allows us to perform detailed diag- 
nostics of the chemical ICM abundances. Consequently, M 87 
is among the best studied objects in the sky. Studies of M 87 
at all wavelengths offer a detailed insight into the balance of 
the ICM cooling and the heating by the central active galactic 
nucleus (AGN). 

The radio map of M 87 shows two lobes of synchrotron 
emission located abo ut 20 - 30 kpc East and South- West of the 
nucleus lOwen et al J 120001) . The AGN and its inner jet prob- 



ably supply the lobes with energy and relativi stic electrons. 
Based on t he previous XMM-Newt on observation Bel sole et alJ 
1200 lb and lMatsushita et al J 120021) found that the ICM in M 87 
has locally a single phase structure, except for the regions as- 
sociated with radio lobes, where the X-ray emission can be 
modelled with two thermal components with temperatures of 
about 2 and 1 keV. The cooler gas fills a small volume com- 
pared to the hotter component and is confined to the radio 
arms, rath er than b eing associated with the potential well of 
the galaxy (Molendi 2002). Based on radio and X-ray observa- 
tions IChurazovetalT i 200 ll) developed a model, according to 
which the bubbles seen in radio buoyantly rise through the hot 
plasma in the cluster and during their rise they uplift cooler X- 
ray emitting gas from the central region out to larger distances. 
This scenario was also studied in numer ical hydrodynamical 
simulations by Briigge n & Raised 120021) who illustrated the 
turbulent mixing of hot, cooler and relativistic plasma initiated 
by the bubble-ICM interaction. This effect probably gives rise 
to the observed two-phase or multi-phase temperature struc- 
ture in the interaction region. Plasma cooling below a tempera- 
ture of 0.8 keV is not observed in the XMM Reflection Grating 
Spectrometers (RGS) X-ray spectra ( Sakel liou et al.l2 002) and 
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neither in the XMM European Photon Imaging Camera (EPIC) 
spe ctra of M 87 dBohringer et all2002l). 

iBohringer etaD l l200ll) and iFinoguenov et al.1 J2002t) 
showed that oxygen has a relatively flat radial abundance dis- 
tribution compared to the steep gradients of silicon, sulfur, 
argon, calcium and iron. The strong abundance increase of 
iron in the core of M 87 compared to that of oxygen indi- 
cates an enhanced contribution of type la supernovae (SN la) 
in the central regions of the galaxy. The flatter oxy gen gra- 
dient was also co nfirmed bv iMat sushita e t alJ J2003I) . On the 
other hand, while iGastaldello & M olendi (2002) confirm the 
strong heavy-element gradient they also find a significant in- 
crease of abundance for light elements, in particular oxygen, 
toward the center of the galaxy which contrad i cts the re- 
sults of |Bohringer et alJ (200l|», IFinoguenov eTafl (|2002) and 
Matsushita et al. M2003l) . 

In this paper we report the results obtained by combining 
two de ep observation with the R eflection Grating Spectrometer 
(RGS, Iden Herder etaD bOOlb aboard XMM-Newton. The 
strongly peaked high X-ray luminosity makes M 87 an ideal 
target for the RGS. The high statistics allows us for the first 
time to reliably determine the abundance of carbon and nitro- 
gen. 

While it is well known that elements from oxygen up 
to the iron group are primarily produced in supernovae, the 
main sources of carbon and nitrogen are still being debated. 
Both elements are believed to originate from a wide range 
of sources including winds of short-lived massive metal rich 
stars, longer-lived low- and intermediat e-mass stars and also an 
early generation of ma s sive stars (e.g.lGus tafsson et al. 1999; 
IChiappini et al . 20031 iMevnet & MaedeiJ 12002b . IShi et alJ 
(2002) find that carbon is enriched by winds of metal-rich 
massive stars at the beginning of the Galactic disk evolution, 
while at the later stage it is produced mainly by low-mass 
stars. They also conclude that nitrogen is pr oduced mainly 
by intermediate-mass stars. iBensbv & Feltzingl l i2006l) find that 
the carbon enrichment in the Galaxy is happening on a time 
scale very similar to that of iron. They conclude that while 
at low metallicities the main carbon contributors are massive 
stars, carbon is later produced mainly by asymptotic giant 
branch (AGB) stars. 

The deep observation also allows us to extract spatially re- 
solved high-resolution spectra and to determine radial temper- 
ature profiles and abundance profiles of several elements out to 
3.5' from the core of M 87. 

For the dist ance to M 87 we adopt a value of 
16.1 Mpc ( To nrv et al.l l2001 ) which implies a linear scale of 
4.7 kpc arcmin -1 . Unless specified otherwise, all errors are at 
the 68% confidence level for one interesting parameter (Ay 2 = 
1). 

2. Observations and data reduction 

M 87 was observed with XMM-Newton on June 19 2000 for 
60 ks and re-observed in January 2005 with an exposure time 
of 109 ks. We process the RGS data from both observations 
with the 6.5.0 version of the XMM-Newton Science Analysis 



System (SAS) following the method described in lTamura et alJ 
(1200 ll) . 

In order to minimise the effect of soft protons in our spectral 
analysis, we extract a light-curve using the events on CCD 9 
of the RGS, outside the central area, with a distance larger 
than 30" from the dispersion axis and cut out all time inter- 
vals where the count rate deviates from the mean by more than 
3<x. After removing the high background periods from the data, 
the net effective exposure time is 40 ks for the first observation 
and 84 ks for the second observation. 

Because the extended emission fills the entire field-of-view 
of the RGS we model the background using the standard RGS 
background model available in SAS. We correct the standard 
response for the RGS instruments by additional, time depen- 
dent large scale effective area corrections found by one of 
us (Jacco Vink; more information can be found on the web- 
site of the XMM-Newton users group 1 ), which are based on 
the standard XMM-Newton calibration sources MRK 421 and 
PKS 2155-304. These corrections will be implemented in the 
forthcoming release of the XMM-Newton SAS software pack- 
age and calibration data files. 

In order to derive spatial temperature and abundance pro- 
files from the RGS we extract spectra from four 1 ' wide regions 
along the cross-dispersion axis of the instrument. We select the 
events from rectangular areas on the CCD strip in the cross- 
dispersion direction (see Fig[3}- 

The observed line emission of the hot gas is broadened by 
the spatial extent of the source along the dispersion direction. 
In order to account for this effect in the spectral modelling we 
convolve the model with the surface brightnes s profile o f the 
galaxy along the dispersion direction jTamura et al.l2004l) . We 
derive the surface brightness profile for each extraction region 
from the EPIC/MOS 1 image in the 0.8 - 1 .4 keV band along the 
dispersion direction of the RGS and convolve it with the RGS 
response in order to produce a predicted line spread function 
(lsf). Because the radial profile for an ion can be different from 
this intensity profile, this method is not ideal. Therefore, we 
introduce two additional parameters, the width and the centroid 
of the lsf which are left free during spectral fitting in order to 
match the observed profiles of the main emission lines. The 
scale parameter s for the width is the ratio of the observed lsf 
width and the width of the intensity profile. 

3. Spectral models 

For the spectral analysis we use the SPEX package 
dKaastra et alJll99r3) . We model the Galactic absorption using 
the hot model of that package, which calculates the transmis- 
sion of a plasma in collisional ionization equilibrium (CIE) 
with cosmic abundances. We mimic the transmission of a neu- 
tral gas by putting its temperature to 0.5 eV. Compared to the 
older s tandard absorption models like lMorrison & McCammonl 
dl983„. the hot model is using continuum opacities from 
IVerner & Yakovl evl d 19951) . Line opacities for most ions are 
from IVerneretalJ d 19961) . To find the best description of the 

1 http://xmm.vilspa.esa.es/external/xmm_user_support/usersgroup/ 
200605 1 8/rgs_calib_eff .pdf 
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emission of the hot plasma in M 87 we use a combination of 
two CIE plasma models (MEKAL) with coupled abundances. 
We also try to model the emission of the hot plasma with 
a differential emission measure (DEM) model with a cut-off 
power-law distribution of emission measures versus tempera- 
ture (wdem). The wdem model appears to be a good empiri- 
cal approximation fo r the spectrum in cooling cores of clus- 
ters of galaxies (e. g. iKaastra et alJl2004l Twerner et al J 120061 
Ide Plaa et alJ I2006T) . The emission measure Y = J n e nndV 
(where n e and «h are the electron and proton densities, V is 
the volume of the s ource) in the wdem m odel is specified in 
Eq. Q adapted from Kaastra "e7aT] ( 12004 : 



dY_ 
~dT 



■ AT i/ a 




elsewhere. 



(1) 



The emission measure distribution has a cut-off at T m i a = 
cT max . For (x — > oo we obtain a fiat emission measure distribu- 
tion. The emission measure weighted mean temperature r mean 
is given by: 



r T ™*dY TdT 

JT min dT 1 U1 

C' m §dT ' 

J L miu ai 



(2) 



By integrating this equation between T m \ a and T max we obtain 
a direct relation between r mea n and T m!lx as a function of a and 

c: 



(1 + l/a)(l -c 1/a+2 ) r 
(2+ l/ar)(l -c l l a+l ) 



(3) 



A comparison of the wdem model with the classic al cooling- 
flow model can be found in Ide Plaa et alJ J2005I) . We note 
that the wdem model contains less cool gas than the classical 
cooli ng-flow model, which is co nsistent with recent observa- 
tions IPeterson et alteOOllEooll) . 

The spectral lines in the MEKAL model are fitted self con- 
sistently. From the fits we obtain the abundances of all elements 
with detected line emission. Throughout the paper we give the 
measured ab undan ces relative to the proto-solar values given 
bv iLoddera (|2003). We note that the values of the fitted ele- 
mental abundances do not depend on the chosen values for the 
solar abundances. We give the values of the elemental abun- 
dances with respect to solar for convenience . The recent solar 
abundance determinations bv lLoddersN2003l) give significantly 
lower abunda nces of oxygen, neon and iron than those mea- 
sured by Ande rs & Grevessd 11989). 



depend on the value for the absorption column density in the 
model. For the best measured va lue of the absorption column 
density N H = 1.8 x 10 20 crrT 2 jLieu et alJI 19961) we obtain a 
power-law with a photon index of T = 2.18 + 0.03 and a 2 - 10 
keV luminosity of (8.6 + 0.5)xl0 40 erg s _1 . 

The AGN was much brighter during the second observation 
than during the first one and also its photon index changed. 
However, the strong pile-up both in EPIC/MOS and EPIC/pn 
does not allow us to determine the properties of the AGN emis- 
sion from the EPIC observations. We have to do it using the 
RGS which is much less susceptible to pile-up. Since the two 
observations were performed with similar position angles the 
thermal cluster emission in the two datasets is the same. We 
determine the parameter values of the cluster emission in the 
first dataset. While we fit the parameter values for the power- 
law emission in the second dataset we keep the parameter val- 
ues of the cluster emission fixed. We find a power-law with a 
photon index of F = 1.95 + 0.03 and a 2 - 10 keV luminosity 
of (1.15 + 0.11) x 10 42 erg s _1 . To confirm this high flux, we 
scale the RGS spectrum obtained during the first observation 
to the exposure time of the second observation and we subtract 
it from the RGS spectrum from the second dataset. This way 
we get rid of most of the cluster emission, which allows us to 
analyze the excess power-law emission from the central AGN 
detected during the second observation. With this exercise we 
confirm our previous results about the high AGN luminosity in 
the second dataset. 

While fitting the thermal ICM/ISM emission we keep the 
power-law parameters fixed at the best determined values. The 
uncertainties in the AGN luminosity introduce systematic un- 
certainties in the emission measure of the thermal components 
and in the absolute abundance values. The uncertainties in the 
photon index of the power-law emission are much smaller than 
those in the luminosity and thus have a smaller influence on 
the best fit parameters of the thermal components. If the lu- 
minosity of the AGN in the second observation is 10% lower 
than the adopted value, then the best fit emission measures of 
the thermal components will be ~ 15% higher and the absolute 
abundance values will be ~ 15% lower. If the luminosity of the 
AGN is 10% higher than the adopted value, then the best fit 
emission measures will be ~ 15% lower and the absolute abun- 
dance values ~ 15% higher. The best fit temperatures and rel- 
ative abundance values are, however, robust and the systematic 
uncertainties arising from the uncertainty in the AGN emission 
are smaller than their statistical errors. 



3.1. The AGN spectrum 

The center of M 87 harbours an active galactic nucleus (AGN). 
Its emission can be well fitted with a power-law. The good de- 
scription of the AGN emission is essential for the determination 
of the absolute abundances in the core of M 87 with the RGS. 
In the first observation we determine the parameter values of 
the power-law by fitting the EPIC spectrum of the core of the 
galaxy. We extract the spectrum from a small circular region 
with a radius of 20" and fit it with a power-law and two ther- 
mal models. The parameter values of the power-law strongly 



4. Results 

4.1. The combined dataset 

Since the difference of the position angles of the two observa- 
tions is 187.3°, the stripes used to extract the spectra from the 
RGS are rotated with respect to each other by only 7.3°. While 
in the first observation the core of the emission falls on the cen- 
ter of the RGS detector, in the second observation the core of 
the emission is offset from the center of the detector by ~1.5'. 
In order to obtain a combined spectrum with the highest possi- 
ble statistics of the core of the galaxy, we extract spectra from 
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Fig. 1. 1st order RGS spectra obtained during two observations 
with different luminosities of the central AGN. The two spectra 
are extracted from a 1.1' wide stripe centred on the core of 
M 87 and are fitted simultaneously. The total useful exposure 
time is 124 ks. The continuous line represents the best fit model 
of the spectrum. 




10 20 

Wavelength (A) 



Fig. 2. Residuals of the fit to the combined RGS spectrum with 
a total exposure time of 124 ks with line emission set to zero in 
the model. We indicate the positions of all the detected spectral 
lines. 

a 1.1' wide stripe centred on the core of M 87 from both obser- 
vations. We fit these spectra simultaneously. The fit to the total 
spectrum is shown in Fig.^ 

We try two different thermal models to fit the emission of 
the hot gas in M 87: a combination of two thermal components 
and a differential emission measure model (see section[3j. The 
results of these two fits are shown in Tabled The reduced ^ 2 s 
of both fits are ~ 1 .4, non of the thermal models can be ruled 
out by our analysis. 

Since the continuum is not very well determined by the 
RGS, absolute abundance values obtained by fitting different 
thermal models are slightly different. However, the relative 



Table 1. Fit results for the RGS spectra extracted from a 1.1' 
wide strip centred on the core of M 87. The ir on abundance is 
given with respect to the proto-solar values bv lLoddeTsI (T2003 ) 
and the other abundances with respect to iron. Emission mea- 
sures (Y = f n e nudV) are given in 10 64 ctrT 3 . 



Parameter 


2T-model 


wdem-model 


Yi 


0.86 ± 0.03 




Y 2 


9.97 ± 0.13 




Y\ydem 




11.56 ± 0.17 


kT { (keV) 


A OA / \ n 1 

0.80 ± 0.01 




kl 2 (KeV J 


1.0/ ± u.uz 




ft J max 




2.38 ± 0.05 


a 




0.47 ± 0.02 


C 




0.22 ± 0.01 


kT 

mean 




1.82 ±0.04 


C/Fe 


0.73 ±0.11 


0.74 ±0.13 


N/Fe 


1.62 ±0.19 


1.62 ±0.21 


O/Fe 


0.60 ± 0.03 


0.59 ± 0.04 


Ne/Fe 


1.40 ±0.11 


1.25 ±0.12 


Mg/Fe 


0.70 ± 0.05 


0.60 ± 0.06 


Fe 


0.98 ± 0.02 


1.06 ±0.03 


x 2 /y 


2737/1905 


2728/1905 



abundances are consistent for both thermal models. In Tabled 
we show our best fit abundances relative to iron, which is, due 
to the strong iron lines in the spectrum, the element with the 
best determined abundance value. The excellent statistics of the 
combined spectrum allows us, for the first time, to determine 
relatively accurate abundance values also for carbon and nitro- 
gen. All the detected spectral lines are shown in Fig.|2] which 
shows the residuals of the RGS spectrum with line emission set 
to zero in the model. 

We note that ISakelliou et all d2002l> use a cooling-flow 
model, with r m ; n = 0.6 keV, in combination with a single tem- 
perature model. By fitting this model we do not obtain a bet- 
ter fit than with the 2-temperature model and the determined 
abundances (also the absolute abundances) are consistent with 
those obtained with t he 2-temperature mod el. The abundance 
values determined by Sakelliou et al are different from 

our values due to a different power-law slope and normalisation 
used to account for the AGN emission. They obtained different 
parameter values for the power-law because they used a differ- 
ent (obtained by a survey at lower spatial resolution) Galactic 
absorption column density of 2.5 x 10 20 ctrT 2 . 

The spectra lack emission from lines that trace gas cooling 
below 0.8 keV (e.g . O VII line em i ssion) , which confirms the 
previous results of ISakelliou et alJ (|2002). The spectra, how- 
ever, allow us to put strong upper limits on the amount of the 
cooling gas. Our upper limit on the emission measure of the 
gas with a temperature of ~ 0.2 keV with the same oxygen 
abundance as the ambient hot gas is Y — 6.9 x 10 61 cm 4 . 
If we assume that all this gas is in a sphere with a radius of 
0.55' then our upper limit of the total mass of the ~ 0.2 keV 
gas is w 10 7 M , which is a; 8% of the total gas in the sphere 
(« 1.2 x 10 8 M o ). 
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4.2. Intrinsic absorption or resonant scattering 

In order to detect intrinsic absorption by the putative cool- 
ing gas in the core of M 87, and/or by the warm-hot inter- 
galactic medium associated with the Virgo cluster (reported by 
iFuiimoto et al we subtract the first observation from the 

second one. This way we get rid of most of the cluster emis- 
sion, which allows us to analyze the excess power-law emission 
from the central AGN detected during the second observation. 
We do not see any absorption lines in the power-law spectrum. 
However, for the O vn absorption column we find an upper limit 
of 10 16 cirT 2 . 

Resonant scattering of line emission is expected to be 
imp ortant in the de nse cores of clusters of galaxies (e.g. 
Gilfanov et allll987l) . Its detection in the RGS spectra of the 
ellip tical galaxy N GC 4636 in the Virgo cluster was reported 



lipti_ 

bv IXu et all J2002I) . The best lines to detect resonant scatter- 
ing in M 87 are the 2p - 3s and the 2p - 3d Fexvn lines at 
17.1 A and 15.01 A respectively. The oscillator strength of the 
2p - 3d transition is substantially higher than that of the 2p-3s 
transitions. Their ratios at different distances from the core of 
the galaxy can reveal resonant scattering of the line at 15.01 A. 
Since these lines originate from the same ion of the same ele- 
ment, the observed spatial dependence of their ratios cannot be 
due to gradients in temperature or elemental abundances in the 
ISM. However, these lines can be well determined only close to 
the core of the galaxy. We determine the ratios of their equiva- 
lent widths in a 30" wide stripe going through the core of the 
galaxy and in the combined spectrum extracted from two 15" 
stripes next to the central extraction region. While the ratio of 
the equivalent widths of the 2p-3s lines to that of the 2p - 3d 
line in the central region is 1.3+0.3, the ratio in the neighbour- 
ing bins is 0.7 + 0.5. This result indicates that 15.01 A photons 
emitted in the core of the galaxy get scattered before they exit 
the interstellar medium. However, the large errors do not allow 
us to make strong conclusions. 

4.3. Radial profiles 

The excellent statistics of the RGS spectra obtained during the 
second observation allows us to extract spectra from four 1' 
wide strips in the cross-dispersion direction of the instrument 
and make radial temperature and abundance profiles up to 3.5' 
from the core of the galaxy. Unfortunately the core of the emis- 
sion is offset from the center of the RGS detector so we could 
not determine the profiles from both sides of the core. 

Since the 2-temperature model describes the data equally 
well as the wdem model and because the 2-temperature model 
may also be well justified by two separate thermal components 
within and outside of the radio lobes in M 87, we fit all the 
spectra with the more simple 2-temperature model. The best fit 
results are shown in Table|5] We note that the extraction regions 
from which we determine the RGS radial profiles are actually 
sampling the South- Western lobe at different distances from the 
core (see Fig- O- 

The radial temperature profiles for the two thermal com- 
ponents are shown in Fig. 0] They both show an increase of 
temperature in the 0.5 - 1 .5' region and they do not show a ra- 




Fig. 3. The RGS extraction regions used for fitting the ra- 
dial profiles overplotted on a temperature map by Aurora 
Simionescu (in preparation). Our radial profiles go along the 
South- Western lobe filled with cooler gas - clearly seen on the 
temperature map. 



> 

0) 



> 

0) 




1 2 
Radius (arcmin) 

Fig. 4. Radial temperature profiles of the two thermal compo- 
nents obtained by fitting spectra extracted from rectangular re- 
gions in the cross-dispersion direction of the RGS. 



dial gradient. The radial distribution of the emission measures 
shows a ~ 3 times larger relative contribution of the cooler gas 
in the central region than in the 0.5 - 2.5' region. 

The radial abundance profiles for oxygen and iron are 
shown in Fig. [5] We see that the radial distributions of oxy- 
gen and iron, the two best determined abundances, are differ- 
ent. Both elements have a peak in the core of the cluster, but 
while the oxygen abundance drops at 1 ' distance from the core 
by ~17% and then its distribution stays flat, the distribution of 
iron shows a gradient all the way out in our investigated region. 

The radial distribution of neon also seems to show a gradi- 
ent throughout the investigated region. But its abundance val- 
ues are more uncertain. In the extraction regions not going 
through the core of the cluster the spatial extent of the source 
along the dispersion direction broadens the lines in the RGS 
and the neon line gets blended with the iron lines. Thus, for 
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Fig. 5. Radial abundance profiles for oxygen and iron obtained 
by fitting spectra extracted from rectangular regions in the 
cross-dispersion direction of the RGS. 

Table 2. Fit results for spatially resolved RGS spectra. We 
model the plasma with two thermal components. Ab undances 
are giv en with respect to the proto-solar values by iLod ders 
(2003). The emission measures are given in 10 64 ctrT 3 . 



Par. 


-0.5/0.5' 


0.5/1.5' 


1.5/2.5' 


2.5/3.5' 


Yi 


0.82 ± 0.04 


0.23 ± 0.03 


0.16 ±0.02 


0.18 ±0.02 


Y 2 


10.83 ±0.17 


8.38 ± 0.08 


6.18 ±0.10 


4.79 ± 0.07 


kT { 


0.79 ± 0.01 


0.81 ± 0.03 


0.69 ± 0.02 


0.69 ± 0.03 


kT 2 


1.69 ±0.03 


1.87 ±0.04 


1.88 ±0.05 


1.66 ±0.05 


C 


0.63 ±0.16 


0.44 ±0.13 


0.30 ±0.16 


< 0.18 


N 


1.64 ±0.24 


0.62 ±0.18 


0.67 ± 0.22 


< 0.23 


O 


0.58 ± 0.03 


0.48 ± 0.02 


0.52 ± 0.03 


0.45 ± 0.03 


Ne 


1.41 ±0.12 


1.17 ±0.13 


0.89 ±0.17 


< 0.23 


Fe 


0.95 ± 0.03 


0.78 ± 0.03 


0.65 ± 0.03 


0.50 ± 0.03 


Scale s 


2.14 ±0.10 


1.10 ±0.05 


1.18 ±0.09 


1.42 ± 0.08 


*7v 


1308/918 


1265/918 


1199/918 


1187/918 



extraction regions further away from the core of the galaxy, the 
systematic uncertainties on the neon abundance get higher. 

The radial distribution of carbon and nitrogen seems to be 
consistent with that of iron, however their values are too uncer- 
tain to draw conclusions. 

We note that the best fit temperature and abundance values 
determined from the RGS are actually the mean values along 
the dispersion direction of the instrument and the true radial 
profiles have stronger gradients than those derived from the 
RGS spectra. 

5. Discussion 

5.1. Chemical enrichment by type la and core collapse 
supernovae 

The ISM/ICM in M 87 was enriched by heavy elements pri- 
marily by core collapse supernovae (SNcc)> type la supernovae 
(SN la) and stellar winds. 



Table 3. Comparison of the observed mass ratios in the core of 
M 87 for five elements to iron with the predicted mass ratios 
assuming that all the elements were produced only by super- 
novae with 86% of all iron produced by SN la. The observed 
to predicted ratio for oxygen is here by definition one. 



Mass ratios 


Obs. 


Pred. 


Obs./Pred. 


M(C)/M(Fe) 


1.32 


0.13 


10 


M(N)/M(Fe) 


0.94 


0.003 


313 


M(0)/M(Fe) 


2.89 


2.89 


= 1 


M(Ne)/M(Fe) 


1.32 


0.36 


3.7 


M(Mg)/M(Fe) 


0.37 


0.19 


1.9 



In order to investigate the relative contribution of SN la and 
SNcc one can use the observed ratios of elements (produced by 
supernovae) with detected line emission to iron to find the best 
fit to the following function: 



/ M(X) \ 
We)j cluster 



f) 



M(X) \ 
M(Fe) J 



(4) 



SN C 



where the expression M(X)/M(Fe) represents the fraction of 
the mass of a given element to the mass of iron as observed in 
the galaxy (cluster), as predicted for the yield of SN la and as 
predicted for the yield of SNcc; / represents the fraction of iron 
produced by SN la relative to the total mass of iron produced 
by the two main types of supernovae. 

For nucleosynthesis products of SNcc we adopt an aver- 
age yield of stars on a mas s range from 10 M Q to 50 M Q cal- 
culated by iTsuiimoto et alJ (1995) assuming a Salpeter initial 
mass function. For nucleosynthesis products of SN la we adopt 
values calculated for the delayed-detonation model WDD2 by 
kwamotoeta i1 lll999l) . 

The best determined are the oxygen and iron abundance. 
Assuming that all iron and oxygen is produced by supernovae, 
from their ratios we estimate that 86% of iron is coming from 
SN la. Then the relative number of SN la contributing to the 
enrichment of the ICM is ~ 40% whil e the relative num- 
ber of SNcc is ~ 60%. We note that iBohringer et all |2005) 
based on the XMM-Ne wton EPIC measurements analysed by 
lMatsushitaetaDd2003l) estimate that about 80% of iron comes 
from SN la. Given the different extraction regions and the gra- 
dient in the iron abundance, this estimate is in a good agree- 
ment with our result. The measured neon abundance is not 
matched well by the supernova models. It is much higher 
than what the models predict (see Table 0, which was pre- 
yiously also observ ed in the cluster of galaxies 2A 0335+096 
(We rner et al.1 2006). The best fit of the magnesium abundance 
value may be influenced by the uncertainties in the effective 
area calibration of the RGS at short wavelengths. Therefore, 
we do not discuss it here. 

The supernovae produce only small amounts of carbon and 
nitrogen. As shown in Table [3] we observe about 300 times 
more nitrogen and 10 times more carbon in M 87 than what 
can be explained by the employed supernova models. The dif- 
ference between the predicted mass produced by all supernovae 
(as calculated from the total mass of iron, relative numbers of 
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Fig. 6. Comparison of the [C/Fe] and [O/Fe] enrichment in 
M 87 with those of the G alactic thin and thick disk stars of 
iBensbv & Feltzinglj2006l) . Thin and thick disc stars are marked 
by open and filled circles respectively. The value for M 87 is in- 
dicated by the asterisk. 



SN la and SNcc and the theoretical yields) and the total ob- 
served mass of these elements can be used to put an upper limit 
on the contribution from stellar winds to the chemical enrich- 
ment of the ICM. 

However, to estimate the total mass of these metals is dif- 
ficult. The long and narrow extraction region of RGS and the 
disturbed nature of the core allow us only to make rough esti- 
mates. If we assume that the measured central abundances are 
typical for the inner 2.5' of the galaxy and w e consider an av- 
erage density of 2.0 x 10~ 2 cirT 3 (based on Matsushita et al. 

than the total mass of carbon and nitrogen in the ISM is 
« 7 x 1 6 M and « 5 x 1 6 M Q respectively. However, the radial 
abundance profiles of nitrogen and carbon suggest abundance 
gradients, which means that since the measured abundance val- 
ues are an average along the dispersion direction of the RGS, 
the true central abundance is probably higher. Because of the 
assumptions, the estimated masses can be considered as lower 
limits. 

5.2. Carbon and nitrogen abundances and enrichment 
byAGB stars 

While it is well known that elements from oxygen up to the iron 
group are primarily produced in supernovae, the main sites of 
the carbon and nitrogen production are still being debated. It 




Fig. 7. Comparison of the [N/Fe] and [N /Ol enrichment in 
M 87 with the Galactic stellar population of Che n et al.l (2000) 
and lShietall J2002I) . The value for M 87 is indicated by the 
asterisk. 



is believed that both elements are being contributed by a wide 
range of sources. The main question is whether their production 
is dominated by the winds of short-lived massive stars or by the 
longer-lived progenitors of AGB stars. 

An analysis of the enrichment histor y of carbon in our 
Galaxy has for example been performed bv lBensbv & Feltzind 
(2006). They find that the ratio C/Fe is fairly constant with a 
decrease of about a factor of 1 .5 over the chemical evolution 
history involving [Fe/H] values from -1 to +0.4 as covered in 
their survey. At the same time the O/Fe ratio decreases by about 
a factor of 5. The general explanation of this large decrease in 
the relative oxygen abundance is an early enrichment by SNcc 
with high O/Fe values and a subsequent primarily iron pro- 
duction by SN la, which go off on much longer time scales 
than the SNcc- In contrast the flat C/Fe ratio then implies that 
the enrichment by Fe through SN la and the C pollution has 
occurred very much in parallel. The similar enrichment time- 
scale suggests that the main sites of carbon production are the 
longer-lived progenitors of AGB stars. 

For our case of the, on average, much older stellar popula- 
tion of M 87 it is now interesting to compare the enrichment 
situation with that of our Galaxy. In Fig. [6] we compare the 
C/Fe and O/Fe enrichment in M 87 with those of the g alac- 
tic thin and thick disk stars of Bensbv & Feltzing (2006). We 
see that while O/Fe is significantly lower in the ISM of M 87 
than in the stellar population of our Galaxy, C/Fe is compara- 
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ble to the lowest C/Fe values detected in the Milky Way. The 
O/Fe ratio in M 87 is so low probably due to the fact that the 
contribution of SN la to the enrichment of the ISM in M 87 
is approximately 3 times as large as in our Galaxy (according 
to Tsuiim oto et al.l l l 1995b the relative number of SN la for our 
Galaxy is Ma/ia+n = 0.13). The low O/Fe ratio also suggests 
that the star formation in M 87 essentially stopped. The C/Fe 
ratio is consistent with the Galactic value which indicates that 
the enrichment by Fe from SN la and by C through AGB stars 
is occurring also in M 87 in parallel. 

In Fig. we compare the N/Fe an d N/O enrich ment with 
a sample of Gal actic stars published bv lChen et alJ J200(j) and 
IShi et alJ J2002I) . The measured N/O ratio in M 87 is higher 
than in the sample of Galactic stars. However, the N/Fe ratio is 
comparable to the highest values detected in this sample. Given 
the old stellar population in M 87 the high N/O ratio in the ISM 
is not surprising. While the oxygen (and probably also some of 
the nitrogen) was supplied by core collapse supernovae early 
in the enrichment history, nitrogen is still being constantly sup- 
plied by intermediate mass AGB stars on a time scale possibly 
similar to the enrichment by iron. 

5.3. Spatial abundance distribution 

While the chemical abundances of all elements have a peak in 
the center of the galaxy, the radial abundance profiles indicate 
that the spatial distributions of oxygen and iron are different. 
The radial distribution of iron, primarily a type la supernova 
product, has a gradient throughout the field of view of the RGS. 
However, the radial distribution of oxygen, primarily a product 
of core collapse supernovae, seems to be flat with a projected 
abundance of ~ 0.5 solar, with only a modest peak in the ex- 
traction region centred on the core of the galaxy. While the flat 
component of the distribution of core collapse supernova prod- 
ucts has been produced by supernovae going off primarily in 
the earlier stages of the evolution of the galaxy, type la super- 
nova still continue to enrich the ISM/ICM in M 87. 

Ou r results confirm the previ ous r esults of lBohring er et alJ 
feOOll) . iFinoguenov et all J2002I) and iMatsushita et all J2003I 
who based on observations with EPIC showed that oxygen has 
a relatively flat radial abundance distribution compared to the 
steep gradients of the heavier elements. While the low spec- 
tral resolution of EPIC combined with the relative closeness 
of the iron lines and calibration problems at low energies left 
some space for doubts about these early results, the RGS with 
its high spectral resolution showed conclusively that the spa- 
tial distributions of oxygen and iron are different. However, the 
systematic uncertainties in the absolute abundance values leave 
still open the possibility of more enhanced abundance values 
in the innermost bin fitted with the RGS, which would mean a 
stronger central peak in the distribution of iron, but also a peak 
in the oxygen distribution. 

6. Conclusions 

We have analysed the temperature structure and the chemical 
abundances in the giant elliptical galaxy M 87 using high- 
resolution spectra obtained during two deep XMM-Newton 



observations with the Reflection Grating Spectrometers. We 
found that: 

- While we confirm the two-temperature structure of the ISM 
we show that a multi-temperature wdem fit also describes 
the data well and gives the same relative abundances as the 
two-temperature fit. 

- The O/Fe and O/N ratios in the ISM of M 87 are lower than 
in the stellar population of our Galaxy, which shows that 
the relative contribution of core collapse supernovae in the 
old stellar population of M 87 to the enrichment of the ISM 
was significantly less than in our Galaxy. 

- The comparison of the C/Fe and N/Fe ratios in the ISM 
of M 87 with those in the stellar population of our Galaxy 
suggests that the enrichment of the ISM by iron through 
Type la supernovae and by carbon and nitrogen is occur- 
ring in parallel and thus the dominant source of carbon and 
nitrogen in M 87 are the intermediate- and low-mass AGB 
stars. 

- From the oxygen to iron abundance ratio we estimate that 
the relative number of core collapse and type la supernovae 
contributing to the enrichment of the ISM in the core of 
M 87 is ~ 60% and ~ 40% respectively. 

- The spatial distributions of oxygen and iron in M 87 are 
different. While the oxygen abundance distribution is flat 
the iron abundance peaks in the core and has a gradient 
throughout the field of view of the instrument, suggesting 
an early enrichment by core-collapse supernovae and a con- 
tinuous contribution of SN la. 
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